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In this talk I present the complete 1-loop perturbative computation of the renormalization constants and 
mixing coefficients of quark and gluon lattice operators of rank two and three whose hadronic elements enter in 
the determination of the first and second moment of Deep Inelastic Scattering Structure Functions, making use 
of the nearest-neighbor improved "clover-leaf" lattice QCD action. 

To perform the huge amount of calculations required for the evaluation of all the relevant Feynman diagrams, 
extensive use of symbolic manipulation languages like Schoonschip and Form has been made. 



1. Introduction 

The computation on the lattice of renormal- 
ization factors is a necessary ingredient to con- 
nect lattice operators and matrix elements to 
their continuum counterparts and extract phys- 
ical quantities from Monte Carlo results. We re- 
port here on the complete computation of the 
renormalization constants and mixing coefficients 
of the quark and gluon operators of rank two 
and of the quark operators of rank three [Q 
that measure the first two moments of Deep 
Inelastic Scattering Structure Functions (SF's). 
We have computed all these constants using the 
Shcikholcslami-Wohlert (SW) 0(a) improved ac- 
tion with csw = 1 in 1-loop perturbation theory. 
Improvement reduces the systematic error 

associated with the finiteness of the lattice spac- 
ing a, which is lowered, for on-shell quantities, 
from O(a) to 0(a/ log a). 

To work out the huge expressions arising for 
each of the numerous Feynman diagrams, it has 
been necessary to use extensively the algebraic 
computer manipulation languages Schoonschip 
and Form. 
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2. Moments of Structure Functions 

The hadronic tensor W^, from which the SF's 
are defined, is given in terms of the hadronic cur- 
rents by 



1 

2^ 
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Using the Wilson OPE near the light-like region, 
A(x)B(0) ~ En^nA* 2 )^ 1 ■ ■■x^OiZ'% N (0), 
the local product of operators can be expressed 
in terms of a set of symmetric and traceless op- 



erators O 



fll—flM 



with vanishing vacuum expecta- 



tion values. The matrix elements of the dominant 
(twist two) operators in the expansion, oj^)..^, 
have the general form 



\p) = A N (/i) p^ ■ ■ ■p flN + traces. (2) 



The formula for the moments of the SF's becomes 



M 

where T x = 2F U F 2 = F 2 /x and T z = F 3 . 
The coefficients Cm are known from continuum 
perturbation theory. From (|j|) we can extract 
(x* -1 ) knowing the corresponding matrix ele- 
ment (p|0^"-/ijv \p)- The latter contains long dis- 
tance (non-perturbative) physics, thus the only 
viable way to compute the moments of SF's is 
with the use of lattice methods. 
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We have considered in our calculations the un- 
polarized SF's, and in particular we have com- 
puted the renormalization constants and mixing 
coefficients of the operators 



0% v = \i)l {ll Du}^ - 
0% v = T. p ^[F{, p F P u } ] 
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0% UT = \i>l {li D v D T} i) — ^ (: 
where is the covariant derivative. 



3. Sheikholeslami-Wohlert improvement 

The matrix elements appearing in (|^) need the 
computation of two- and three-point correlation 
functions || . To reduce discretization errors the 
SW improved action is used [Q. Adding to the 
Wilson action the "clover-leaf" SW term 
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and performing on the spinor fields in the opera- 
tors 0^„ and 0^„ T the "rotations" 
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leads, in on-shell matrix elements |3j], to a can- 
cellation of all terms that in the continuum limit 
(<7o ~ 1/loga) are effectively of order "a" QJ^]. 

The systematic error related to the lattice dis- 
cretization drops in this way from order a to order 
a j log a. Numerically this results in a substantial 
improvement. In fact, while the magnitude of the 
order a terms is about 20-30 %, the magnitude 
of the order a/ log a terms turns out to be about 
5-10 %. 

4. Renormalization constants 

The renormalization constants connect the 
bare lattice operators, O(a), to finite operators, 
0(/i), renormalized at a scale y,: 



0\y) = Z lk ^a)O k {a). 



(5) 



In the flavor Singlet case there is a mixing be- 
tween quark and gluon operators of the same rank 



that have the same conserved quantum numbers. 
We then write: 



0« = Z qq Oi + Z qg 03 



(6) 



03 = Z gv Oi + Z gg 03, 
and in this case all elements of the mixing matrix 



(q\0«\q) (g,a\0«\g,cT) 
(q\0<>\q) (g,a\Oa\g,a) 



(7) 



have to be computed. 

To this mixing (already present in continuum 
QCD) the lattice regularization adds additional 
mixings, induced by the breaking of (Euclidean) 
Lorentz invariance 0. It is possible, by a care- 
ful choice of the Lorentz indices, to simplify at 
least partially the lattice mixing pattern. How- 
ever, the higher the moment the more compli- 
cated the mixing pattern becomes fjl|8| . In partic- 
ular, the non-Singlet operators 0' 12 j and 0? 12 3} 
are multiplicatively renormalizable on the lattice 
in the quenched approximation, but another use- 
ful second moment operator, Oq IS = ^?4m — 

~l(°{i22} + °{433})> is 110t - 

The breaking of Lorentz invariance forced us 
to develop special computer routines to correctly 
perform the Dirac algebra on the lattice. They 
play a key role in our codes which are designed 
to automatically carry out all the steps of the 
necessary algebraic manipulations, starting from 
the elementary building blocks of each Feynman 
diagram. 

5. Some results 

Putting Z — 1 + g% /167r 2 (7 log ay+B), one gets 
for the -B's and j's of the first moment operators 
the results of Tables 1 and 2. The numbers 
reported here for the B's are for r = 1 and differ 
slightly from the ones published in the first paper 
of jlj. They have been re-computed by us and 
cross-checked against the works in || . 

Simulations have been performed in the past 
with the unimproved Wilson action [|]J|] . Within 
errors, Monte Carlo results are consistent with 
experiment. No simulation has been as yet per- 
formed in the improved case. Values of the renor- 
malization constants in the improved theory are 
somewhat larger than in the standard Wilson 
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Wilson 


Improved 


(b£) w 


-3.165 


B f 


-15.816 


0.208 




-0.743 


-5.817 


-4.044 


-2.168 


-6.084 




B aa = 


-15.585 





Table 1 

The lattice constants B for the first moment; B? g 
corresponds to the quark loop contribution to the 
gluon propagator. 
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B qq 


vertex 
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ANOMALOUS 




self-energy 


-1 


DIMENSIONS 




total 


- 34/9 




16/3 


B qg 


- 4/9 


Iqg 


4/3 


B gq 


- 22/9 


Igq 


16/3 


Bgg 


- 4/3 


Igg 





Bi 
"gg 


- 10/9 


igq 


4/3 



Table 2 

The continuum constants B in the MS scheme 
and the anomalous dimensions for the first mo- 
ment. 



case. As an example, we give here a selection 
of the results for quark operators (for the com- 
plete results see [Q). At j3 = 6.0 and for r = 1 
one gets 



(Of 



12}) 



1.027 (O? ^ w 



{0\ l2} f = 1.134 (0\ 12} f 

(Of 123} ) W = 1-160 (Of 123} ) w 

(dl 23} f = 1.252 (0\ 123} Y 

(Odis) W = I [1-184 (O a ) w + 1.156 (0 B ) W ] 

(OhsY = I [1-331 (Oa) 1 + 1-187 (0 B ) 1 } , 

where the first of each pair shows the Wilson re- 
sults (for which we are in agreement with ||) and 
the second the improved results. 

Finally, we want to mention that results now 
exist also for the polarized SF's, though limited 
to the Wilson case B. 
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